Temperature dependance of the tunneling density of states in sub-micron planar 

metal / oxide / graphene junctions 
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We present tunneling measurements of sub-micron metal/insulator/graphene planar tunnel junc- 
tions up to room temperature. We observe a gate independent gap, as previously observed only by 
low temperature STM[Y. Zhang et al., Nat. Phys. 4, 627 (2008)]. No gap appears at temperatures 
above l^QK, which is four times smaller than the theoretically expected Tc, from the accepted mean 
field model[T. O. Wehling et al. Phys. Rev. Lett. 101, 216803 (2008)]. We show that taking into 
account an additional vibrational effect of out-of-plane phonon soft modes the gap may disappear 
from the measurements at temperatures much lower than the calculated Tc. 



Graphene [l| is famous for its linear dispersing elec- 
tron and hole bands which touch at two nonequivalent 
points K and K' at the corner of the Brillouin zone, so 
it is gapless. The Fermi level of graphene lies at these 
Dirac points. Hence, it was a surprise, that instead of 
a linear V-shape conductance characteristic, a ~ lOmeV 
gap, in the tunneling spectra in STM experiments was 
observed 0, [sj at low temperatures. The characteris- 
tics were symmetric and gate independent gap around 
zero bias voltage, which means it opens up at the Fermi 
level rather than at the Dirac point. The gap was inter- 
preted as due to an electron-phonon interaction that 
couples Dirac electrons at the K point (tt band) with 
nearly free electrons at the F point (cr band), where the 
gap opens up. 

We present tunneling measurements of metal / Al- 
oxide / graphene sub-micron planar junctions up to room 
temperature. The broad temperature range is made pos- 
sible by the planar geometry. Our data presents the same 
tunneling gap as seen so far only using STM at low tem- 
peratures. Unexpectedly we find that it disappears at 
^ 150K, a temperature 4 times smaller than expected 
from the accepted model [1]. We suggest an addition to 
this model, whereby we take into account the vibrational 
effect of out of plane soft mode phonons. This may cause 
the gap to disappear at a temperature much lower than 
the mean field Tc- 

We fabricated sub- micron metal / insulator / graphene 
junctions, as is shown in Fig. [21 The flake is outlined. 
We used commercially available single layer flakes [5] ex- 
foliated on highly doped Si substrate with a 300nm ox- 
ide. The flakes have a typical mobility of a few thousand 
crn^/Vs. The patterns were defined using standard e- 
beam lithography. In the first step metallic Ti/Au con- 
tacts were evaporated using e-gun. In the second step 
we evaporated 20A of Al in O2 gas at partial pressure of 
IQ-'^Xorr with the substrate cooled to 77K. The oxide 
layer covers the whole region of the flake. In the final 
step we evaporate an additional electrode, on top of the 




-0.1 0.1 

V[Volt] 



'Electronic address: ' hacohe@post.tau.ac.il| 



FIG. 1: Differential conductance as function of temperature. 
The data is shifted up for clarity. The intersection of the two 
orange lines, drawn on samples B's 5K curve, exemplify the 
extracted gap value. 



oxide and across the flake, this serves as the top elec- 
trode of the junction. We fabricated junctions with Pb 
electrodes (sample A) and with Ti/Au electrodes (sample 
B) and obtained similar results with both. We performed 
DC differential conductance measurements using Keith- 
ley 6221-K2182A DeltaMode. 

The low temperature differential conductance spectra 
of the junctions, as shown in Fig. [U shows a gap of 
~ 7QmeV. The gap is constant for all gate voltages, 
measured from -50V to 50V. The gap is characterized by 
a relatively constant conductance at low bias voltage fol- 
lowed by a rapid rise in the conductance. We chose the 
bias voltage of this intersection point to represent the 
gap value, as illustrated by the two orange lines, drawn 
on sample B's 5K curve in Fig. [TJ Fig. [2] shows the ex- 
tracted gap closes at ^ 150K. 

The low temperature gap is the same, in size, as ob- 
served by STM [1,11 , which shows that it can be probed 
using sub-micron planar junctions. The experimental Q 
and theoretical 0| interpretation is that this is a phonon 
assisted tunneling gap, where the phonon is an intrin- 
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FIG. 2: a) Extracted gap value as a function of temperature, 
as described in the text. Green and blue symbols are the 
values extracted from the data. The red squares are extracted 
from the simulated data. Fig. [S] Inset - Calculation of the 
order parameter as a function of temperature, showing that 
the calculated gap closes at 590K. The gap value is b) 
Scanning Electron Microscope image of sample A. The red 
dashed lines outline the flake. 



sic out of plane mode with energy of QlmeV which hy- 
bridizes the TT (Dirac electrons) and a (nearly free elec- 
trons) bands. As a consequence a gap in the a band 
(nearly free electrons) opens up and a rigid shift of the 
energies in the tt band (Dirac electrons) occurs. 

To calculate the temperature dependence of the gap 
we take the model interaction Hamiltonian given by T. 
Wehhng et. al. [1] 



Hm = A^(4+gC, + cl^qdq){ak + alfc). (1) 

q.k 

Where dq is the annihilation operator of an electron 
in the a band, in this band the electrons are nearly 
free electrons and may be approximated by a flat band 



with energy eq = 3.3ey. 



is the annihilation oper- 



ator of an electron in the tt band(Dirac electrons), au 
is the annihilation operator of an out of plane phonon 
with momentum k and A is the electron-phonon cou- 
pling strength. We apply a BCS-like mean field ap- 
proach and define an order parameter, /S.^ = ^k,k' < 
cl.,dk' >■ The self energies of the problem are then 
E± = ± 
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The excitation spectrum has a gap of size 

nearly free electron band, and a rigid shift of - 
energies of the Dirac electrons, in exact agreement with 
the results of T. Wehling et. al. [3]. We can now proceed 
and write the gap equation: 



Afe — ^ Vk,k' 

k' 



Afe' 



-.{fd{E^) - UE+)) (2) 



Where fd is the Fermi function. If we take Ak.k' = A 
and Vk,k' — V, such that the gap is momentum inde- 
pendent, we can numerically solve the equation, where 

we used the experimental gap value — — 67meV. 

The calculated temperature dependance of the gap ^ is 
shown in the inset of Fig. [2l The gap closes at 590K. In 
the calculation wc took a momentum independent inter- 
action, yet taking some form of a momentum dependent 
Vk^k' will correspond to stronger coupling and only raise 
the calculated Tc. In contradiction to the calculated Tc 
of 590K, the observed gap disappears from the measure- 
ments at about 150K. Thermal smearing will come into 
effect only at temperatures above 400K and cannot ex- 
plain the disappearance of the gap at about 150K. 

To see how this difference can be reconciled we must 
first understand why, when tunneling to the tt band, the 
measurement is sensitive to a gap in the cr band. Fol- 
lowing T. Wehling et. al. Q, the tunneling conduc- 
tance is proportional to the density of states of each 
band multiplied by the density of the electronic wave- 
function at the position of the counter electrode dl/dV ~ 
I'Srl'^NriE) + |^'i^|2Ar^(£'). On the graphene sample 
VPrp is negligible compared with I^'k^. But l^'ifp de- 
cays away from the surface, while l^'rP increases up to 
a distance of ~ lA from the surface and beyond that ex- 
ponentially decays into the vacuum. The two densities 
decay on different length scales, and for distances larger 

than lA their ratio is j.^'-'L = e^'^'', where h is the dis- 
tance in angstroms between the graphene sheet and the 
electrode. Accordingly when the electrode is far away 
from the graphene sheet the tunneling conductance is 
mostly into the a band, but close to the sample, it is into 
the TT band. At the typical distances used in STM and 
planar junction experiments, tunneling is mostly into the 
cr band, which is why a gap is observed. 

The idea behind our suggested explanation is that at 
high enough temperature the excited out of plane phonon 
soft modes will give rise to vibrations large enough to ef- 
fectively bring the graphene sheet close enough to the 
electrode, such that tunneling will be directly to the tt 
band and the gap will not appear in the data, despite 
it's existence. The amplitude A{u!) of a phonon of fre- 
quency u) is ^ pT^j where p is the density and V is the 

volume of the sample. Soft mode phonons, at low en- 
ergies, have a high density of states. Also the existence 
of such modes was considered 0-01 as the mechanism 
behind the additional unexpected rise in resistivity in 
graphene, above 150K on 5*102 @, and above lOK in sus- 
pended graphene Q . This suggests that the contribution 
of these low energy excitations cannot be overlooked. 

M. Mohr et al. @ showed that graphene has an out 
of plane soft mode band (ZA). To evaluate the effect of 
the vibrations we calculate y/ < A'^{u!) > as a function of 
temperature. 
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FIG. 3: a) The amplitude of the wavefunction of the Dirac 
electrons at the position of the tunnel electrode plotted as 
function of temperature, b) Simulated data realizing the tran- 
sition from a a band density of states to a tt band density of 
states. Simulated data is equally spaced in temperature from 
20K to 180K and shifted up for clarity. The lowest temper- 
ature (20K) simulated plot corresponds to the function we 
used to represent the g band and the highest temperature 
(180K) simulated plot corresponds to the function we used to 
represent the tt band. 



< A^{uj) >= J A^{Lj)fb{oj,T)D{oj)duj. (3) 

Where fb is the Bose function and D{uj) is the phonon 
density of states. The density of graphene is p — 
7.6 10-'^^. We took the disper sion relation as k(uj) = 
O.lSw'''^®, by fitting the out of plane acoustic phonons 

band (ZA)[9|, where ui is in meV and fc is in 4. The fit 

A 

is suitable for energies up to SOmeV, beyond these ener- 
gies /ti(cLi, T)D{ljj) is practically zero. Since the integrand 
diverges at zero energy we must use a cutoff parameter 

^cutoff- 

We now plug the contribution of the vibrations into the 
conductance. The sum of the densities |^r|^ + l^ifl^ = 1 
and the ratio, away from the surface, is p^J^ = e^^''^ so 
the fraction of the tt band electrons is, 



' ' 1 + gl.7(/io-\/<A2>) ^ ' 

Where is the distance in angstroms between the 
electrode and the graphene sheet at zero temperature. 
The tunneling probability to the tt band will substan- 
tially increase as the vibrations approach the value of 
(2nm). The dependence of \^k\^ on the temperature, 
shown in panel 'a' of Fig. [3l shows there is a defined 
transition temperature, where the contribution of l^rl'^ 



is diminished and l^'xp becomes dominant. In this cal- 
culation we used ujcutoff = 2.2 x lO^^meV . Since the 

ratio |^^|2 = e^'^'* holds only away from the surface, 
and closer to the surface the ratio decays rapidly to zero, 
the transition beyond lAOK will be more rapid in temper- 
ature than illustrated in panel 'a' of Fig|3l We extracted 
the gap from the simulated data, the results are plotted 
in Fig. El 

Since the oxide thickness {ho = 2nm) is experimentally 
known, the only parameter in the calculation is the cutoff 
frequency uj cutoff, which was chosen to give a transition 
temperature similar to the one observed experimentally. 
It corresponds to a cutoff wavelength Xcutoff — 2Ap,m, 
which is a reasonable number slightly smaller than the 
dimensions of the physical flake. The actual flake size 
doesn't define \cutof / since all the measurements show a 
similar transition temperature. Our calculation is based 
on the graphene modes Q as measured in graphite. Yet in 
our experiment the flake is partially pinned to the 5*^02 
substrate and has the Al-oxide on top. While in both 
cases (graphite and "oxide sandwich") the graphene is 
weakly coupled to the top and bottom layers and prob- 
ably has a very similar phonon spectrum, the different 
environments may alter the spectrum of the vibrations. 
More importantly a close inspection of the SEM image. 
Fig. [21 shows a non-uniform and possibly porous struc- 
ture on the scale of tens of nanometers which belongs to 
the Al-oxide layer. Such a non-uniformity means that 
at a fixed temperature the vibrations will be different at 
different positions, so in the thinner regions of the Al- 
oxide layer the transition will occur when the amplitude 
of the vibrations are substantially smaller than /iq and 
that the thicker regions may contribute less to the tran- 
sition. Therefore the transition temperature is actually 
a result of some averaging over the junction region. In 
a STM measurement, where there is no oxide separating 
the graphene and top electrode, so the transition tem- 
perature may be lower. 

In summary we have shown the ability to fabricate 
planar tunnel junctions to graphene that are sensitive 
to the tunneling density of states. So far this has only 
been achieved using STM. Our technique opens up the 
possibility to tunnel into graphene in sub-micron sized 
devices. We measured the phonon assisted tunneling gap 
over a broad temperature range. Unexpectedly we found 
that the gap disappears at '^ISOK, although the Tc cal- 
culated from the accepted mean field model is 4 times 
larger. We suggest that the early disappearance of the 
gap may be due to the effect of out of plane soft mode 
phonons, whose thermal excitations bring the wavefunc- 
tion of the TT band electrons, which have no gap, close to 
the counter electrode. 
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